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Abstract 

Baryons with a single heavy quark are being studied experimentally at present. 
Baryons with two units of heavy flavor will be abundantly produced not only at future 
colliders, but also at existing facilities. In this paper we study the production via 
heavy quark fragmentation of baryons containing two heavy quarks at the Tevatron, 
the LHC, HERA, and the NLC. The production rate is woefully small at HERA and 
at the NLC, but significant at pp and pp machines. We present distributions in various 
kinematical variables in addition to the integrated cross sections at hadron colliders. 



The constituent quark model has been remarkably successful in describing the ob- 
served hadronic states. This includes mesons with one heavy and one light quark (i.e., B 
and D mesons), baryons with one heavy and two light quarks, and the J/ip and T mesons, 
which contain a heavy quark and heavy antiquark. Recent experiments have made great 
progress in the observations of the c and b baryons |2| and their decays ||, [§. The con- 
stituent quark model predicts the existence of baryons containing two heavy quarks (cc, be, 
or bb) and one light quark, and we are on the verge of obtaining the necessary experimen- 
tal sensitivity to observe these states. The energies necessary to produce these particles 
are already reached; the difficulty remaining is in their reconstruction. These states have 
in general a large number of decay modes, so that their observation and a measurement 
of their properties will require a large number of them to be produced. This difficulty is 
increased by the fact that the production rates at e + e~ colliders are extremely small, so that 
the identification of these particles must take place in the messier environment of hadronic 
collisions. 

Were the heavy quark of infinite mass, the two heavy quarks would be bound in a 
point-like diquark, and the light degrees of freedom would "see" this diquark as a static 
antitriplet color source. In this limit, the interactions of the light degrees of freedom with 
the heavy diquark are quite similar to interactions of the light degrees of freedom with the 
b quark in a B meson[|5], |J. For realistic heavy quark masses, this simple picture is not yet 
valid. In particular, the hyperfine splittings of these baryons are not yet well described by the 
heavy- mass limit 0. The states may nonetheless be accurately treated as a combination of 
light degrees of freedom and heavy although not-pointlike diquark. The interactions between 
the two heavy quarks are analogous to the QQ system familiar from J / ip and T spectroscopy. 
At short distances, the interaction will be dominated by one-gluon-exchange, although with 
a color factor 2/3 from the QQ interaction rather than 4/3 from QQ, and at long distances 
it will be confining 0. 

The production of these states is reliably calculable as a hard process. One takes 
the full set of Feynman diagrams for the production of two pairs of heavy quarks and heavy 
antiquarks and requires that the momenta of the two heavy quarks be nearly equal, then 
projects out the correct quantum numbers for the relevant Q1Q2 state. This is similar to the 



1 



calculation of the hard production of heavy quarkonium states|9|], though more complicated 
as an additional QQ pair is needed. On the other hand, because the quarks are heavy, frag- 
mentation functions for a heavy quark to produce a doubly heavy diquark can be predicted 
in perturbative QCD|TO], |TT|. The fragmentation description of the Q1Q2 production process 
is not valid where the masses of the heavy quarks become important, that is, for small p T 
or small energies, but for high-energy colliders such as LEP or Tevatron it should be appli- 
cable. The calculation is similar to that described in [O for J/ip production. The form of 



the matrix element is different in the two cases, requiring a somewhat changed method of 
projection onto the proper spin states. For the production of a QQ bound state, the standard 
procedure is to replace the heavy quark production matrix element uOv = Tr (Ovu) with 

the trace Tr (OPssJ, where Pss z is the proper projection onto the spin of the QQ state |I2 



This procedure is, however, not simply applicable to the Q1Q2 bound states due to the form 
of the matrix element, which contains not u and v as in the production of QQ but u\ and 
«2- The projection procedure described above is nothing more than a clever method of per- 
forming the sum over the quark spins. This sum can also be performed by summing over 
helicity amplitudes with the proper Clebsch-Gordan coefficients. Although the procedure is 
somewhat different, the calculations are identical in their essentials, and the fragmentation 
function is simply proportional to that for a QQ pair. The differences reflect the changed 
color and statistical factors. The fragmentation function for Q\ to produce a spin-1 diquark 

(Q1Q2) is 

Afc->«ww(*.M = K + 2m *)) = ™ 12lR{QlQ i m2 a 2 s (2rn 2 )F(z), (1) 

F ( z ) = 77^ 2 ~ 2 ^( 3 - 2r ) + 3 ^( 3 - 2r + 4r " 2 ) 

w 4(1 - z + rzf L v ' y ' 

-2^ 3 (1 - r)(4 - r + 2r 2 ) - z 4 (l - r) 2 (3 - 2r + 2r 2 )] , (2) 

where r = m 2 / (m 1 +m 2 ). The factor iV 12 = (1)2 for (un)equal quark masses, and reflects the 
presence of identical fermions in the final state. For production of a spin-0 state, (Q1Q2)', 
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the function F(z) must be replaced by 



F'(z) 



rz(l — z 



■2 



6 - 18z(l - 2r) + z 2 (21 - 74r + 68r 2 



12(1 -z + r^) 6 

-2^ 3 (1 - r)(6 - 19r + 18r 2 ) + 3z 4 (l - r) 2 (l - 2r + 2r 2 )l . (3) 



As in the case of hard quarkonium production, the requirement of nearly equal momenta 
and the projection onto the proper Q1Q2 quantum numbers will modify the p T spectrum of 
the diquark; it will fall more rapidly than that of single heavy quark production. For heavy 
quarkonium production, the faster fall-off in p T of the hard production process means that 
fragmentation production will dominate for sufficiently large p T despite the suppression of 
the fragmentation production by powers of a s [13j. For the Q1Q2 states, no hard produc- 
tion processes of lower order in a s exist. We assume that the doubly heavy diquarks will 
always hadronize into baryons with two heavy quarks, so that the fragmentation function 

-DQi^baryon = I?Q 1 _ f (Q 1 Q 2 ). 

A Q1Q2 diquark in the ground state has a symmetric spatial wavefunction. Because 
the Q1Q2 must combine with a light quark to produce a colorless baryon, the diquark must 
be a color antitriplet (antisymmetric) state. For identical quarks the wavefunction must 
satisfy Fermi-Dirac statistics, and thus the cc and bb ground state diquarks have spin 1. The 
be states are not restricted by statistics, and we denote the spin singlet (triplet) state by be' 
(be). The spin triplet can hadronize into a spin-1/2 or -3/2 baryon, 3^ and S£ c , respectively; 
the spin singlet can produce only the spin-1/2 "E' bc . The probabilities for the fragmentation of 
heavy quarks into these baryons are: c — > H cc , S* c , about 2 x 10~ 5 ; b — > E' bc , about 4 x 10~ 5 ; 
b — > E bc ,El c , about 5 x 10~ 5 . The remaining probabilities are suppressed by ~ (m c /m,b) 3 
and therefore approximately two orders of magnitude smaller . 

In addition to these probabilities, it is desirable to have predictions of the p T (or other 
kinematical variables) distributions for the production of these states. The fragmentation 
functions given above at low momentum scale /io must be evolved to higher scales fi in order 
to be of use in calculations for high energy colliders: 

da(A + (Q1Q2) +X) = Yl£ dz da ( A + i{ rf> v) + X, /1 2 ) D^ {QlQ2) (z, fi 2 ). (4) 



The evolution occurs via an Altarelli-Parisi equation: 

^ 2 -Q^ D i^(QiQ2)(z,fi 2 ) = 7^J2J z — P ij( z /V) D j^(QiQ2)(y^ 2 )- (5) 

which does not change the overall probabilities given above. In the evolution of the frag- 
mentation functions via (|5]), we include only the Pqq splitting function. The fragmentation 
of gluons to heavy diquarks is suppressed relative to that of heavy quarks by a s , and may 
be neglected. A separate numerical evolution of the fragmentation functions to each z and 
H needed in the course of a numerical simulation would be an extremely computer-intensive 
approach; instead we generate values of the fragmentation function for a reasonably-spaced 
set of \i and z values and interpolate to the desired points. The unevolved fragmentation 
functions depend on |i?(0)| 2 , the non-relativistic radial wavefunction at the origin, the heavy 
quark masses, and a s (fj,o). These input parameters are given in Table |I[ 




Table 1: Parameters used to calculate the unevolved fragmentation functions 



The J/ip and T decays to lepton pairs give the wavefunctions |i?(0)| 2 for these 
states relatively accurately |16| and independently of phenomenological potential models. 
The QCD-corrected expression for the leptonic decay of heavy quarkonium states 



2„2 



r(v 



\Rv(o)f 



! 16 , > 



(6) 



is used to extract the T and J /if} radial wavefunction at the origin. There is no information 
on the B c meson to determine the wavefunction for the be states, and the annihilation of this 
state into lepton pairs would not exist in any case to provide the measurement. However, 
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there is an empirical relationship among leptonic decays of the 0, J /if) and T mesons, r(V — > 
e + e~) ~ 12eg keV, which indicates that |i?(0)| 2 //i 2 is approximately constant. 

The one-gluon exchange contribution for Q1Q2 differs from that for Q±Q 2 by a relative 
color factor 1/2. In a Coulomb potential, therefore, the relation \Rj/^(Q)\ 2 = 8|i? cc (0)| 2 must 
hold. The heavy quarks are, however, not in a pure Coulomb potential, so that we must 
rely on phenomenological potentials to better estimate the radial wavefunction. While many 
functional forms for the potential provide quite good fits to the heavy quarkonium data, the 
fitted potentials all have very similar shapes in the region of r ~ 1 GeV -1 . We use 1/2 
the Q1Q2 potential for the Q1Q2 potential, and find the wavefunction at the origin using a 
numerical solution of the Schrodinger equation. We compare in Table [| the results of this 
calculation for two familiar potentials with the values extracted from the J /if) and T decay 
rates. 





Indiana fl4|| 


Richardson |15| 


data|22j 


cc 
cc 
be 
be 

bb 
bb 


(0.95 GeV) 3 
(0.67 GeV) 3 
(1.18 GeV) 3 
(0.82 GeV) 3 
(1.72 GeV) 3 
(1.15 GeV) 3 


(0.93 GeV) 3 
(0.65 GeV) 3 
(1.18 GeV) 3 
(0.81 GeV) 3 
(1.88 GeV) 3 
(1.23 GeV) 3 


(0.98 GeV) 3 
(1.96 GeV) 3 



Table 2: Radial wavefunctions at the origin, |i?(0)| 2 , for QiQ 2 and Q1Q2 states for two 
interquark potentials compared to values extracted from data, see (6). 

Doubly heavy baryons will be produced at all present and future accelerators where 
there is sufficient energy. The results presented here include both baryon and antibaryon 
production. We consider first the e + e~ colliders LEP, LEPII, and NLC. At LEP, the only 
subprocess for the production of a heavy quark Q is e + e~ — ► Z, with the Z decaying to QQ. 
The production rates in this case are quite small, a few events per year. For LEPII and NLC 
the additional processes e + e~ — > W + W~ where one W decays to a c quark and e + e~ — > ZZ 
where one Z decays to cc or bb also contribute (depending on the LEPII energy, of course). 
Despite the additional subprocesses, of which W + W~ dominates for NLC, the event rates 
remain woefully small. At LEPII, only a few events/year can be expected, and a 60 fb _1 /yr, 
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500 GeV linear collider will produce only of order 30 events per year. The production rates 
at 77 colliders and cy colliders were also calculated, including backscattered and Weizsacker- 
Williams photons and the contributions from resolved-photon processes. The cross sections 
for production of these particles at these colliders are also hopelessly small. 

The results at HERA are similar. The main subprocesses here are 75 — > QQ and 
iQ ~~* 9Q- Assuming the design luminosity of 200 pb _1 /yr, we again find that the event 
rate is rather small, of order 30 events/year. The resolved photon processes are likewise 
negligible. 

The situation is considerably more hopeful at hadron colliders. The relevant subpro- 
cesses are gg — ► QQ, qq — > QQ, gQ — > gQ, and qQ — > qQ (in the latter two subprocesses, 
q and Q stand for both quarks and antiquarks), although the main contribution is from the 
gluon-fusion subprocess. For the Tevatron, with an integrated luminosity of 100 pb _1 /yr, 
we anticipate of order 8 x 10 4 events/year {p T > 5 GeV, \r}\ < 0.5), while the LHC, with 
100 fb _1 /yr, should produce of order 1.3 x 10 8 events per year(p T > 10 GeV, \r]\ < 0.5). Total 
production cross sections per unit rapidity for the various diquarks at the Tevatron and the 
LHC (operating at both 10 TeV and 14 TeV) are given in Table |3|. The rates for production 
of the frc-type states are similar to those expected for the B c meson, where rates of 10 4 _B c 's 
have been predicted at the Tevatron (p T > 10 GeV) and 10 7 at LHC (p T > 20 GeV)|iT|. On 
the other hand, the production rate for cc-type states is significantly smaller than that for 
J /if} production at the Tevatron which is expected to be 5 x 10 6 (p T > 5 GeV) [p~3"l . 





Tevatron 


LHC (10 TeV) 


LHC (14 TeV) 




430 pb 


330 pb 


470 pb 




145 pb 


220 pb 


330 pb 


^bc + ^bc 


215 pb 


350 pb 


490 pb 


^bb -r ^bb 


16 pb 


27 pb 


36 pb 



Table 3: Production cross sections per unit rapidity in the central region. 

Some results are given in Figs, [i] and |2|. Fig. [I] shows da / dp T / di]\ v=Q vs. p T . A 
kinematical cut of 5 GeV for the baryon p T has been imposed here. Fig. |^ gives the rapid- 
ity distribution of Q1Q2 baryons, again with p T > 5 GeV. The rapidity distributions are 
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rather flat in the central region, so that da/dp T may be estimated from da/dp T /dr]\ ri= o by 
multiplying with the ^-region desired. 




Pt [GeV] 

Figure 1: p T distributions (at small rapidity) 
for doubly heavy baryon production at the 
Tevatron. 
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Figure 2: Pseudorapidity distributions for 
doubly heavy baryon production at the Teva- 
tron. 
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Figure 3: p T distributions (at small rapidity) 
for doubly heavy baryon production at the 
LHC, 14 TeV. 
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Figure 4: p T distributions (at small rapidity) 
for doubly heavy baryon production at the 
LHC, 10 TeV. 



Figs. |3| and £| show da/dp T /dr]\ ri= o vs. p T for the LHC operating at 14 TeV (Fig. ^|) 
and 10 TeV (Fig. f|). A kinematical cut of p T > 10 GeV is placed on the Q1Q2 baryon. The 
production rate at the lower center-of-mass energy is slightly lower, but even at the lower 
energy, the LHC will produce Q1Q2 baryons copiously. The rapidity distributions are again 
relatively flat in the central region (\rj\ < 3), and event rates for any specific detector can be 
estimated by simply multiplying the results in Table |3] by the r\ coverage of the detector. 
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In order to estimate the uncertainties due to the choice of factorization and fragmen- 
tation scales, we vary these scales (all chosen to be equal) by a factor of 2 from our nominal 
choice fx = p T . The effect of varying the scale on da / dp T / dr]\ v=0 vs. p T can be seen in Fig. ||, 
and the effect on da/dr] vs. 77 can be seen in Fig. |[ The processes shown are representative 
of all the processes. The behavior seen in Fig. ^ is somewhat counter-intuitive in that one 
expects the rate with /1 = p T /2 to be larger than that for \x = p T . However, the fragmenta- 
tion functions are cut off at low fif rag - This causes the drop at low p T for this choice seen 
in Fig. [| and the lower cross sections in Fig. || We use MRSA(T^] parton distribution func- 
tions. The effect of changing the parton distribution set choice is negligible — far smaller 
than the effect of varying the choice of \x. Other sources of theoretical uncertainty include 
the parameters used for the unevolved fragmentation functions, primarily |i?(0)| 2 where dif- 
ferent phenomenological potentials give about 10% differences in the result, unknown QCD 
corrections to the parton-level cross sections, the fragmentation function and the evolution 
equations, and unknown relativistic corrections to the initial fragmentation functions. 




20 40 60 80 100 



Pt [GeV] 

Figure 5: p T distributions (at small rapidity) 
for bb and be production at the Tevatron for 
different choices of scale \i. 
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Figure 6: Pseudorapidity distributions for bb 
and be production at the Tevatron for differ- 
ent choices of scale \i. 



Recently there has been some concern that the fragmentation approximation may 
not work as well as anticipated for B c production. This assertion is based on a comparison 
between the full calculation and the fragmentation approximation for B c in 77 collisions 1 19]. 
There is, of course, a significant difference between 77 and gg initial states, and the con- 
volution of the gluon distributions with the parton-level subprocess can have a large effect. 



S 



The calculation of gg — *> B C X is much more complex than that of 77 — > B C X, but is nec- 



essary to understand the validity of the fragmentation approximation. The authors of 



went on to study B c production at hadron colliders [IT]. They found that the fragmentation 



approximation still differs from the full calculation at the parton level (cr[gg — > bb{-^ B C X)] 
vs. o\gg — > B c bc]) except at high p T , but after the convolution of a with the gluon distribu- 
tion functions, the full calculation and the fragmentation approximation agree quite well for 



p T > 10 GeV both at Tevatron and LHC energies. Fig. 6 of JT7| shows that the agreement 

between the full calculation and the fragmentation approximation is acceptable down to 

pmm _ 5 Q e y a t the Tevatron and p™ m = 10 GeV at LHC. As the fragmentation functions 

and the production mechanisms are extremely similar for B c and heavy diquark production, 

our calculation of the production of doubly heavy baryons should also be accurate in this p T 
range. 



The decay modes of doubly charmed baryons have been studied in |20| using SU(3) 
flavor symmetry. This approach relates the decays of these particles, and could be used to find 
relations among the decays of the bb baryons and among the bc baryons. The transitions 
between the doubly heavy baryons || and some decays of the bc baryons have been 
discussed in the heavy quark limit. The lifetimes of c and b quarks are not very different 
(t a ~ 0.2 x 10~ 12 s while r A ~ 1.07 x 10~ 12 sp2]), so that the weak decays of both b and c 



are important when considering the decays of the bc baryons. bc baryons will decay weakly 
(both hadronically and semi-leptonically) to single b quark baryons as well as doubly charmed 
baryons. Both of these possibilities have been discussed in the literature^, ^3], |24|. Either a 
single b baryon or a doubly charmed baryon traced back to a displaced vertex using a vertex 
detector will be a clear sign of doubly heavy baryon production |E5|. In the case of bb baryons, 
the production of a doubly charmed baryon (possibly with same sign di-lepton!), with both 
6-quark decays tagged by a vertex detector, would be an impressive signal. Unfortunately, 
the bb baryons are produced less frequently than the other doubly heavy baryons, and two 
semi-leptonic decays will severely reduce the event rate. 

In conclusion, the production via fragmentation of baryons containing two heavy 
quarks has been calculated in the fragmentation approximation. The fragmentation approx- 
imation reproduces well the full calculation of B c production at hadron colliders for the 
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p T studied here ||17|| ; it is expected that the fragmentation approximation will work well for 
doubly heavy diquark production as well. The production rates of these baryons at e + e~ 
colliders and at the ep collider HERA are found to be negligible. The situation is much 
better at hadron colliders, with approximately 8 x 10 4 events/yr expected at Tevatron and 
1.3 x 10 8 events/year at LHC In addition to predictions for total production cross sections 
(\r}\ < 0.5 and p T > 5(10) GeV at the Tevatron (LHC)), the distributions da/dp T jdr\ 1^=0 
and da/drj for p T > 5(10) GeV at the Tevatron (LHC) are studied. The detection of these 
particles and measurements of their properties will provide an experimental challenge due 
to the large number of their decay modes, but will be a rich testing ground, e.g., for HQET 
and QCD potential models. 
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